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ABSTRACT 


This study analyses the trend for initial detection times using both passive and 
active sonar during submarine-on-submarine operations. Specifically, it simulates a 
nuclear powered submarine (SSN) searching for a diesel submarine in an environment 
where the SSN has a speed advantage and active sonar detection ranges exceed passive 
sonar detection ranges. The simulation uses a mover-sensor discrete event application of 
SIMKIT, developed by Professor Arnold Buss. 

The simulation results show that initial detection times of a search follow an 
exponential trend as a function of SSN speed, diesel submarine speed, detection ranges, 
ping interval and detection probability. As a result, as detection ranges continue to 
decrease due to increases in sound quieting technology, initial detection times during a 
submarine search will increase exponentially. This can render a passive sonar 
prosecution ineffective when combating a modern diesel submarine. Should an SSN use 
active sonar, initial detection times can be significantly reduced, especially if combined 
with an effective search path. The threat to the SSN of using active sonar can be 
mitigated by judicious consideration of ping interval and search speed with detection 
probability and active detection ranges. All values used to arrive at the conclusions 
stated are notional, and no classified information sources were consulted as part of this 


work. 
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DISCLAIMER 


The reader is cautioned that computer programs developed in this research may 
not have been exercised for all cases of interest. While every effort has been made, 
within the time available, to ensure that the programs are free of computational errors, 
they cannot be considered validated. Any application of these programs without 


additional verification is at the risk of the planner. 
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EXECUTIVE SUMMARY 


There has been an increase in production worldwide of modern diesel submarines. 
Nations have found that these ships can provide a relatively inexpensive ability to both 
defend local coastlines and provide a credible anti-surface warfare threat. Traditionally, 
the nuclear submarine fleet (SSNs) has been trained to passively detect and track 
opposing submarine forces. As these passive detection ranges continue to decrease with 
modern sound-silencing techniques, the ability to detect these diesel submarines 
diminishes, and it must be determined how to properly utilize SSNs to search for these 


threats. 


This thesis uses discrete event simulation to observe trends in both passive and 
active sonar detection. The simulation recreates a submarine—on-submarine operation, 
with an SSN searching for a diesel submarine until initial detection occurs. The 
simulation is performed using SIMKIT, a software simulation tool developed by 


Professor Arnold Buss of the Naval Postgraduate School. 


The initial analysis focuses on performing a random search of a diesel submarine 
using first passive sonar detection, then active sonar detection. This analysis shows that 
initial detection time follows an exponential relationship between detection ranges and 
submarines’ operating speed. Additionally, ping interval and probability of a ping 
resulting in detection contribute to the exponential effect when using active sonar. Due 
to this relationship, it can be concluded that as passive detection ranges continue to 
decrease, the time needed to locate a target submarine will increase exponentially. It also 
shows that as the detection range of active sonar exceeds that of passive sonar by a factor 
determined by the conditions in which active sonar is implemented, it rapidly reduces the 


time to initial detection. 


The second part of the analysis explores how the search path affects the time to 
initial detection when using active sonar. It was found that a "lawnmower" search path-- 
that is, boundary to boundary sweeps as the submarine travels down an area--did not 


improve performance above that of a random search. Though diesel submarine speed is 


XV 


relatively slow when submerged, it has sufficient speed in some instances to cross the 
perimeter of an area previously searched by the searching submarine, especially with 
long ping intervals and lower detection probabilities. However, a "spiral-out" search path 
was found to significantly reduce initial detection time by making it less likely for a 
diesel submarine to enter previously searched areas. This is due to the perimeter of the 


searched area being initially small and increasing linearly as the search progresses. 


The final analysis involves how ping interval is affected by other ships’ 
parameters during a search. If the Commanding Officer deems it necessary to minimize 
the ping interval to be able to minimize the risk to his own ship through the use of active 
sonar, other factors will change the effects of the exponential contribution of ping interval 
to detection time. With a long detection range, a reduction in the SSN's search speed as 
well as a moderate to high detection probability will allow for a lengthening of the time 


between pings with a relatively small increase in detection time. 
The conclusions of the thesis are: 


e The initial time to detect a target is an exponential function based on 
detection range, speed of the participants, and in the case of active sonar, 
ping interval and detection probability. As a result, as the passive 
detection range of diesel submarines continues to decrease, the average 
detection time will increase exponentially. A point exists where a passive 
prosecution is no longer feasible, and either other assets need to be 
employed or an active search could be performed to maintain a reasonable 


detection time. 


e The speed of a submerged diesel submarine is limited; however, it is 
sufficient to make many search patterns no more effective than a random 
search if the diesel can enter previously searched areas. A "spiral-out" 
path can provide an additional decrease in detection time, as it provides an 
initially small perimeter that the diesel can cross in order to enter this 


previously searched area. The simulations suggest that an active search 
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can yield detection times at a fraction of passive searches, but real world 


exercises are needed to determine the breakpoint at which this occurs. 


If the Commanding Officer deems it necessary to reduce ping interval to 
minimize the threat of being counter-detected, the following factors (listed 
in order of priority) will mitigate the effect of the exponential increase in 


detection time: 


1. A long active detection range 
2s Traveling at slower search speeds 
3: A moderate to high detection probability 
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I. INTRODUCTION 


This chapter summarizes the thesis. It includes the problem statement, the focus 


of the research performed, and the expected outcome of the analysis. 
A. PROBLEM STATEMENT 


When a submarine is conducting search operations in an attempt to detect another 
submarine, it is generally conducted using passive sonar. This method detects the target 
submarine once its radiated noise level exceeds that of the ambient noise of the 
environment by way of a receiver mounted on the searching vessel. Another method that 
is available is the use of active sonar, which emits a signal into the environment that 
reflects off a target submarine and is received by the searching submarine. This places 
the searching submarine at risk, as this signal is also generally heard by the target, who 
may then evade or act against the searcher. With modern diesel submarines, however, 
the radiated noise it emits is generally very small, making passive detection difficult. 
During events in which detection of these submarines is vital or the threat of the diesel 
submarine is relatively low, active sonar is a tactical option that must be considered 


during the prosecution of diesel submarines. 


Currently, many nations are developing or purchasing modern diesel-electric 
submarines that are extremely quiet when submerged and operating on the battery. 
Recent technological advances allow them to remain submerged for extended periods of 
time at low speeds. Their low radiated noise has reduced the passive detection ranges to 
the point that the effectiveness of passive prosecution is substantially reduced, especially 


in areas of high ambient noise. 


The current fleet of U.S. nuclear powered attack submarines (SSN) is equipped 
with active sonar. As stated previously, active sonar effectively broadcasts the location 
of the SSN, placing it at risk if the diesel submarine is able to obtain the bearing of the 
SSN based on the radiated noise. The slow speed of the diesel submarine makes evasion 


difficult, but, if equipped with a capable weapons system, can allow the diesel to fire 


upon the SSN. Given this risk, active search is still a viable option for specific high- 
threat situations involving a diesel submarine or against those diesels that are less capable 


of successfully engaging a modern nuclear-powered submarine. 


B. RESEARCH FOCUS 


This thesis uses discrete event simulation to analyze the use of passive and active 
sonar to search for and detect diesel submarines by SSNs. The proliferation of diesel 
submarines throughout the world poses significant challenges to the U.S. dominance of 
the sea. Though unable to travel and operate in vast expanses of the ocean as effectively 
as nuclear powered submarines, diesel submarines are adept at patrolling local littoral 
waters and provide a credible Anti-Submarine Warfare (ASW) and Anti-Surface Warfare 
(ASUW) threat. Though incurring risk, an active sonar prosecution is viable in several 
scenarios to reduce the time to locate diesel submarines. Examples of these include but 


are not limited to: 


e The deployment of significant numbers of diesel submarines by a nation in 
which, though a state of war does not exist, tensions exist. This assumes 
the opposing nation has not authorized the diesel submarines to engage 


any tracking submarines. 


e Diesel submarines are equipped with weapons that cannot effectively 
engage SSNs. For example, their torpedoes are not effective against SSN 
countermeasures or are designed/designated for a high priority surface 


target. 


e A nation sponsoring terrorism has supplied a diesel submarine to terrorists 


which intend to conduct operations in the immediate future. 


The lack of use of active sonar in recent years has likely led to atrophy in 
Commanding Officers’ experience in effectively using active sonar, specifically in the 
knowledge of the benefit (reduction in detection times) that can be gained, and the proper 


employment of active sonar. 


An analysis of the trends of detection times using both passive and active sonar 
provides a guide to the Commanding Officer such that future options have a greater 
chance of success in detecting diesel submarines. The analysis also provides information 
on how to reduce the risk of counter-detection by minimizing the frequency at which a 


ping needs to be emitted using an active sonar search. 
cy OUTLINE 


This thesis is organized into five chapters. Chapter I is divided into four sections. 
The first section describes the problem statement. The second section explains the focal 
point of the research. The third section gives a brief overview of each chapter in this 


thesis. The fourth section explains the expected research outcomes. 


Chapter II explores the different participants that are considered in this study. A 
description of the current U.S. SSN fleet is given. Selected diesel submarines common in 


the world and their capabilities are also discussed. 


Chapter III defines the discrete event simulation and scenario. In the first section, 
the simulation program is explained, along with the parameters and the assumptions 
designated for the simulation. The second section discusses the scenario used as the basis 


of the specific simulation. 


Chapter IV covers the analysis of the simulation results, and contains four 
sections. The first section describes the passive sonar search of a diesel submarine. The 
result of the first section is used as the baseline for further evaluation of following 
searches. The second section analyzes the active sonar search and compares it to those of 
the passive sonar search. The third section analyzes a simulation of an active sonar search 
using different search paths. The fourth section examines ping interval, and how altering 


ping interval and associated search parameters affects initial detection. 


Chapter V summarizes the results obtained from simulations. The conclusions 


provide a guideline for future sonar searches, and areas of further research are identified. 


The Appendix includes results from the simulation runs and the regressions 


performed on the data. 


D. EXPECTED OUTCOME OF RESEARCH 


The thesis explains the trend of initial detection rates as a function of the 
parameters of passive and active sonar searches. The results can be used as a guide in 
future sonar searches by Commanding Officers to provide a better understanding of how 
search type, parameters and actions will affect the time to detect a target submarine. It 
can also be used by senior naval personnel to determine which assets can be effectively 


utilized to conduct a search for diesel submarines. 


Hl. RELEVANT SUBMARINE DESCRIPTIONS 


In order to establish context, this chapter provides a brief introduction of the 
submarines considered in the analysis. This chapter is not intended to thoroughly 
examine each submarine and its capabilities, but rather to explain relative performances 


and provide a background for the model. 


A. SUBMARINES OF THE UNITED STATES 


1 Bs Los Angeles Class Submarine 


The Los Angeles Class submarine is the mainstay of the U.S. Navy's attack 
submarine fleet, with 46 currently commissioned.!2 This will remain in effect in the 
foreseeable future, as its replacement, the Virginia Class submarine, is currently at a 


production rate of two per year.3 
2. Seawolf Class Submarine 


There are three Seawolf Class submarines in the U.S. Navy, with no further 
construction planned. This submarine class was designed to conduct Cold-War era 
submarine warfare, and is faster, better armed, and has an improved sonar system 


compared to the Los Angeles Class submarine.* 


! Navy Fact File “United States Navy Fact File Attack Submarines — SSN”, 
http://www.navy.mil/navydata/fact_display .asp?cid=4100&tid =100&ct=4 (accessed Feb 2, 2008). 


2 Stephen Saunders, Commodore, RN, editor, Janes Fighting Ships 2007-2008, (Surrey, UK: Jane’s 
Information Group Limited, 2005), 884. 


3 DefenseLink News Transcript "Briefing on the Virginia Class Submarine Contract", 
http://www.defenselink.mil/transcripts/ transcript.aspx ?transcriptid=3224 (accessed Feb 2, 2008). 


4 Saunders, 883. 


3. Virginia Class Submarine 


Virginia Class submarines will eventually replace Los Angeles Class submarines 
as the primary attack submarine in the U.S. fleet. Though not as fast and well-armed as 
the Seawolf Class, its modular design allows it to act in multiple roles in most ocean 


environments.° 


B. DIESEL SUBMARINES 


1. Kilo Submarine 


Countries that maintain a Kilo class submarine in service include China and Iran. 
This dual-hulled submarine has an at-sea endurance of 45 days, and can perform multiple 
missions, from mine-laying to anti-naval operations. It is one of the significant assets 


sold by Russia to foreign nations.® 
ze Song Submarine 


In October 2006, a Song class submarine surfaced within five miles of the aircraft 
carrier Kitty Hawk (CV-63). The Song Class is the first modern-era Chinese-built 
submarine, and is armed with wake homing torpedoes and anti-ship cruise missiles. The 


incident emphasized the threat that quiet diesel submarines pose to carrier strike groups.’ 





5 Saunders, Janes Fighting Ships 2007-2008, 882. 


1 


6 GlobalSecurity.org, "Kilo Class Submarine - People's Liberation Army Navy, 
http://www. globalsecurity.org/military/world/china/kilo.htm (accessed Feb 2, 2008). 


7 GlobalSecurity.org, "Song Class - People's Liberation Army Navy," 
http://www. globalsecurity.org/military/world/china/song.htm (accessed Feb 2, 2008). 
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HI. SIMULATION MODEL AND SCENARIO 


This chapter explains the model that is used in the thesis. The program used in 
the simulation, the construction of each entity, how parameters are implemented in the 
model, and the scenario are discussed. Assumptions that are made in the model are 


explicitly stated. 


A. SIMULATION MODEL 


1. SIMKIT 


The simulations performed in this thesis were developed using the Java-based 
program SIMKIT, developed by Professor Arnold Buss of the Naval Postgraduate 
School. The discrete event simulations execute an event graph, in which subsequent 
actions are appended onto an event list. These actions are continuously ordered and 
performed chronologically, and time periods without events do not consume 
computational overhead. The benefit of discrete event simulation is the ability to perform 
multiple simulations in a fraction of the time of equivalent time-step simulations, since 
the model calculations are performed after each event rather than after each time 
increment. In this thesis, over 36,000 distinct simulations were performed, representing 


millions of total simulated hours. 


2. Submarine Representation 


The submarines are represented as uniform linear movers. A more complicated 
representation is not necessary in this model for two reasons: First, the time scale of this 
model (based on the size of the simulation area and speed of the vessels involved) do not 
require a detailed analysis of acceleration, deceleration, advance, and transfer. Second, 
the ability to change depth is not analyzed in the model because any change in depth can 
be represented as a new two-dimensional problem, as depth changes should only have an 


effect on the sonar detection range and the probability of receiving a signal. 


A uniform linear mover starts its move at position x at time ¢, and starts to move 
with velocity v. As a result, its location at time f is x+(f- f,)v. As a result, the position of 
the mover does not have to be explicitly stored at all times in the model, but rather can be 


calculated as needed. 


A separate mover manager maintains the scheduling of all movements in the 
simulation associated with the specific mover. As such, this manager is able to store a 


set of waypoints that the mover is to follow during the conduct of the simulation. 


Though submerged diesel submarines can operate above 10 knots, this places a 
significant drain on the battery and cannot be sustained for long periods of time. 
Therefore, the simulation assumes the diesel submarine travels between two and five 
knots. The SSN is assumed to have a speed advantage over the diesel submarine, and 
will travel between five and fifteen knots (fifteen knots giving a 3:1 speed advantage over 


the diesel's highest speed). 


For the random passive and active sonar searches, each submarine is randomly 
generated in the simulation area. The Kilo performs random maneuvers during all 
simulations. This occurs even if the Kilo is pinged by active sonar. Though the Kilo 
may determine that it is being pinged, it is assumed that it is unable to resolve the bearing 
from which the ping originates. The SSN also maneuvers randomly during the initial 
active searches to provide a baseline for further analysis, in which definitive search paths 


during active sonar prosecution are performed. 


a Sonar Detection Representation 


Both active and passive sonar systems rely on the principle of signal excess. In 
all underwater detection, detection range relies on multiple factors depending on the 
target position, searcher position, and ocean characteristics. Bottom depths, ocean 
bottom contour, composition of the ocean floor, target aspect, and several other factors all 


have an effect on detection range. 


In passive systems, a receiver is mounted to the hull of the searching submarine. 


As radiated noise leaves the target submarine, its intensity decreases as a result of 
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spreading and attenuation losses. A receiver mounted on the searching submarine 
receives the signal and is processed to be evaluated by a sonar operator. If the radiated 
noise level is sufficiently greater than the ambient noise level, detection will occur. This 
detection relies on the operator being able to isolate and recognize the radiated noise as 
that of the target submarine. In passive sonar detection, continuous monitoring of the 
passive sonar system leads to very quick detection once the source (diesel) travels within 
the range that the signal-to-noise ratio is above that required for detection. As a result, 
passive sonar can be reasonably represented as a cookie-cutter sensor. The detection 
ranges vary widely between the types of submarines and operating environments, and the 
simulations assume a notional detection range between 0.5 and 3 nautical miles (nm). 
With active sonar, the SSN will emit a ping (radiated noise) that will travel, reflect off a 
target, and return to the SSN. If the return signal is sufficiently strong, it will be read as a 
detection event. Pings also suffer from spreading and attenuation losses. Though this 
limits the range that active sonar can be effectively used, the source signal generated is 


such that its effective range will be greater than passive range. 


Clutter such as sea mounts and biologics can distort, reduce, or confuse the signal 
such that detection is not certain even if the target is within the range required for 
detection. Increased training can assist in the ability to discriminate a return signal as the 
target, but cannot eliminate all errors. As a result, when a ping occurs and the target is 
within the detection range of the SSN, there is an associated probability that the ping both 
returns to the SSN and is then recognized as the target. It is assumed that the 
Commanding Officer would recognize and not operate active sonar in a poor acoustic 
environment. It is also assumed that active sonar detection ranges will be greater than 
passive detection ranges. The simulation uses an active detection range between 5nm and 
15nm, a ping interval between 0.01 and 1 hour, and a detection probability of 0.3 to 1.0. 
Baffle zones are not considered in the simulation, in order to maintain this thesis’s 
unclassified status. This should not affect results, as the speed advantage of the SSN 


makes detection of the diesel submarine very unlikely to occur in such an area. 


In this simulation, the sensor is a separate object that is attached to the mover. 
The sensor maintains a list of contacts that have been detected, as well as maintaining the 
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parameters needed to calculate if a detection event occurs. A separate referee is 
implemented that will determine whether a detection occurs as a mover enters the range 
of the sensor. Similarly to a mover manager, the referee calculates and stores the times at 
which a mover will enter detection range, be detected, remain undetected, and exit 
detection range. As the movers change course, these times are recalculated and updated 


on the event list. 
4. Parameter Input 


As submarine speeds, detection ranges, ping intervals, and detection probabilities 
all act as a continuous range, there are virtually unlimited combinations of these 
parameters. In order to efficiently ensure that the possible range of these combinations is 
considered, a Nearly Orthogonal Latin Hypercube (NOLH) was used to input the 
parameters of the simulation. This provides uniform samples for the marginal 
distribution of each single input, with the main effects being nearly orthogonal to each 
other (allowing for greater distinction of these effects).8 Each data point in the NOLH 


design was replicated for 200 trials (simulations). 


A general representation of the model is shown in Figures 1 and 2. These were 
obtained directly from the passive detection model to aid in the visual representation of 
the model. The blue square represents the SSN, and the red square represents the Kilo. 
The yellow circle surrounding the SSN represents the sensor's range--in this case either 
passive or active sonar. Each submarine will continue to maneuver until the Kilo 
penetrates the yellow detection range as shown in Figure 2. In the case of passive sonar, 
detection will then occur. If the sensor represents active sonar, detection will only occur 
when a ping is scheduled and a drawn random number is then compared to detection 
probability. An example of the coding used for this simulation is included in the 


Appendix. 


8 Thomas M.Cioppa,"Efficient Nearly Orthogonal and Space-Filling Experimental Designs for High- 
Dimensional Complex Models" (master’s thesis, Naval Postgraduate School, 2002) . 
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Figure 1. Visual Representation of Initial Model Generation 
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Figure 2. Visual Representation of Initial Model Detection 


3. Model Intent 


The discrete event simulation construction is optimistic when compared to real 
world scenarios. Any best possible outcome for a flexible model that can cover multiple 
current and future scenarios requires having an active detection range that is significantly 
high. As a result, the area that the simulation covers is a 125nm x 125nm area. It should 
also be noted that the initial detection times are not predictive, but rather show trends and 


relationships between factors and tactics. 
B. SCENARIO DESCRIPTION 


A U.S. Los Angeles class submarine has been tasked with locating and tracking a 
Kilo submarine prior to the commencement of a carrier strike group exercise. The Kilo 


has recently been deployed and is believed to be located in an area in which a carrier task 
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group will be conducting exercises. The exercise will be conducted in a 125nm x 125nm 
box, and will not include ASW operations. There has been an increase in tension with 
the nation that has deployed the Kilo, and due to the sensitive nature of the exercise it is 
desired that any external presence be identified prior to the decision to initiate the 
exercise. The Kilo has orders to remain in the area unless detected due to the high 
priority its government places on the mission. The Kilo is unable to resolve the bearing 
of any active sonar it detects, and as such will move randomly in the area to avoid 


developing any pattern that can be exploited by searchers. 
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IV. SIMULATION ANALYSIS 


This chapter contains the results and analysis of the various simulations 
performed in the thesis. The initial simulations include baseline passive and active sonar 
searches that provide a baseline for further analysis. Later simulations analyze the 
performance of search patters in relation to detection times. The final portion of this 
chapter examines how to effectively utilize active sonar when a minimum ping interval is 


desired. 
A. RANDOM PASSIVE SEARCH 


The simulation involves both the SSN and the diesel submarine moving randomly 
in the designated area. Though a random search path is not optimal, it has historically 
been shown to be a good representation of reality.2 This is due to the fact that target 
speed, though small, is sufficient to provide some evasion capability. This causes most 
exhaustive search paths to perform no better than a random search. Both submarines are 
generated at a random position in the simulation area. The SSN travels at a constant 
speed between 5 and 15 knots, the Kilo travels at a constant speed between 2 and 5 knots, 
and passive detection range varies from 0.5 to 3 nm. The Kilo does not respond if 
pinged, which may be interpreted as an inability to localize the direction of the ping 
received. Each data point was replicated 200 times, and the results of the simulation are 


shown in Table 1. 


9 Alan Washburn, "Models Based on Detection Rate" Search and Detection, 4th ed., (Institute for 
Operations Research and the Management Sciences, 2002), 2-5. 
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(Knots) (knots) (nm) Detection (hours) (hours) 

| a Sic 348.916 344.70302 
| = 487.147 | __ 455.51209 | 
| Bile 1627.281 | __ 1652.0335 | 
| ni 812.296 
| | 374.088 
| is{ 410.684 
| 1 | 248.216 
| = 303.318 
c—— = 452.372 
| __ 689.644 
| 4 | 707.931 
| 14 | 208.467 
a BE 320.392 
| al 622.731 
| | 745.226 
| la | 1763.783 
| ol =e i 808.765 
Po C—~idEC Overall Mean 643.015[ 
PCE Overall Std Dev 596.585] 


Table 1. Random Active Search Results 


The mean value of detection times over the range of submarine speeds and 
detection ranges is 643 hours. The standard deviations of the results are not constant, 


suggesting heteroscedasticity. Additionally, the values of the standard deviations are 


range, also suggests that detection time looks similar to general exponential decay. 
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Figure 3. Passive Random Search Detection Time 


Since the detection time is not a linear function, a linear regression is not 
appropriate for this data. As a result, an intrinsically linear regression is performed by 
taking the natural log of initial detection times. The results of the regression are shown in 
Figure 4. The initial model that includes all factors shows the p-value for the Kilo speed 
is 0.148. Though common sense would dictate that the Kilo’s speed will affect the time it 
will be detected, the model suggests that the other factors have a much greater effect, to 
the point of making the Kilo's speed negligible. A subsequent regression, shown in Figure 
5, was performed without the Kilo’s speed. Removing this factor did not change the 
amount of variation that is explained by the regression (R squared remains at .2324), 
suggesting that the Kilo’s speed did not explain any more of the variation in detection 
times. The F statistic p value is essentially zero, implying a statistically significant 
model. Both passive detection range and SSN speed have individual t statistic p-values 
of essentially zero, which strongly suggests that these factors are important in predicting 
initial detection times, with passive detection range being the most important factor. The 
R squared states that the model only explains 23% of the variability in detection times, 
but this was expected due to the random submarine location generation and random 
courses of the submarines involved over such a large area. The formula resulting from 
the regression is: 


ly 


Initial Detection Time(hrs)= 


2916.1] ey a8 - 0.1*(SSN Speed, knts) - 0.567*(Passive Detection Range, nm) 


The conclusion drawn from this equation is that as detection range decreases, the 
time to initial detection will increase exponentially. This can be partially offset by 
increasing SSN search speed, but it cannot be overcome. As a result, as the detection 
range of diesel submarines continue to decrease, it is likely that a passive search using an 
SSN will be ineffective and will not find the diesel submarine in any reasonable period of 


time. 
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Observations (or Sum Vgts) 3400 
| Analysis of Variance 

Sum of F-significance 
Source DF Squares Mean Square F Ratio 
Model 2 971.2218 323.741 342.8128 p value 
Error 3396 3207.06S8 0.944 Prob>F 
Cc. Total 3399 4173.2876 <.0001* 
| Parameter Estimates 
Term Estimate StdError tRatio Prob>jt| 
Intercept 7.9734544 0.085331 93.50 0.0000* 
SSN Spd (hours) -0.100753 0.005499 -13.32 <.0001* 
Kilo Spd (hours) -0.02S5637 0.017448 -1.47 0.1418 
Pass Rng (nm) -0.567381 0.022402 -25.33 <.0001* 
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| Sorted Parameter Estimates | 


Term Estimate StdError t Ratio Prob>|t| 
Pass Rng (nm) -0.5673381 0.022402 -25.33 <.0001* 
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Figure 4. Full Regression of Passive Random Search 
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B. RANDOM ACTIVE SEARCH 


This SSN will now use active sonar to search for the Kilo in the area. The Kilo 
and SSN will move randomly in the area, with the SSN ping interval and detection 
probability being set at a constant value. The simulation stops upon initial detection of 
the Kilo. The detection range of active sonar is 5 to 15 nm as compared to the 0.5 to 3 
nm of passive sonar. The lower bound was selected because active detection range will 
always be greater than passive detection range. This will take into account the 
Commanding Officer’s use of active sonar only if its range is significantly greater than 
passive detection range. The higher bound was selected to give a 5:1 advantage of 
maximum passive detection range and maintain the area of simulation at a reasonable 


size. SSN and Kilo speed and limitations will remain the same as in the passive search. 


The results, shown in Table 2, again suggest an exponential trend of average 
initial detection time. The mean time to initial detection is 110.74 hours, as opposed to 
643 hours for the passive random search. Figure 6 also suggests that an exponential trend 


is reasonable. 
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Kilo Active | Pinginerat Average Initial Fea 
Speed Range Detection Ping Interval Detect Time Standard 
(knots) (nm) Probability (hours) (hours) Deviation 

200.81185 

71.03755 

133.14265 

35.69445 

67.2388 

78.49365 

95.82585 

28.6538 

72.67396985 

87.25525 

50.9184 

36.0078 

501.518 

63.06215 

75.19385 

65.37125 

219.14805 

110.743 | 


Overall Std Dev 


Table 2. Random Active Search Results 
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Figure 6. Active Random Search Detection Times 


A full regression performed using the natural log of initial detection time, as 
shown in Figure 7, shows an F statistic p value of essentially zero, and all factors have a t 


22 


statistic p-value less than 0.05. The R squared states the model explains only 20% of the 
variability in the simulation. Though yielding 3% less than the passive search, more 
factors are involved, and there is some probability the Kilo will not be detected even if it 
is pinged inside of the active detection range. The most important factor is active 
detection range, followed by ping interval, detection probability, SSN speed, and Kilo 
speed, respectively. The regression results in the following equation to predict initial 


detection time: 


Intial Detection Time(hrs)= 
A374" *(Active Range, nm)+1.02*(Ping Interval, hrs)-0.86*(Detect Prob) - 0.055*(SSN Spd, knts) + 


0.047*(Kilo Spd, knts) 
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| Analysis of Variance 

Sum of 

Source DF Squares Mean Square F Ratio 
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Error 3394 5990.0675 1.765 Prob>F 
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Estimate StdError tRatio Prob>jt| 


Term 

Intercept 6.078588 0.163146 37.26 <.0001* 
SSN Speed (knots) 0.05499 0.007589 -7.25 <.0001* 
Kilo Speed (knots) 0.0474707 0.023716 2.00 0.0454* 
Active Range (nm) 0.173232 0.007447 -23.26 <.0001* 
Det Prob 0.860985 0.106181 8.11 <.0001* 
Ping Intw 1.0165373 0.075458 13.47 <.0001* 
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Full Regression of Active Random Search 


This regression suggests that the exponential effect of ping interval, though 
important to the model, can be delayed (up to a point) without a significant increase in 
initial detection time. This is due to the diminishing returns inherent in the exponential 
distribution. As ping interval is shortened, it will yield sequentially less reduction in 


initial detection time. This is further explored in Section D. 
cC. SEARCH PATH EFFECTS ON TIME TO INITIAL DETECTION 


The previous model was performed with the SSN conducting a random search. 
This is not an optimal search configuration, and any Commanding Officer performing a 
search would undoubtedly choose a specific patrol path to conduct it. There are multiple 
methods in which a sonar search can be conducted. This thesis considers two such search 


paths. The paths were selected based on capabilities of the SSN and ease of use. 


The lawnmower path of the SSN begins the search in the southwest corner of the 
area and conducts a west-—to-east sweep, traveling some distance north, and then 
conducting an east to west sweep. This pattern is continued until the entire area is 


searched. 


The next path approximates a "spiral-out" search. The SSN will start at the center 
of the area, and then travel outward in areas of expanding squares until the entire area is 


searched. 


A "spiral-in" search path is also possible but was not included in the 
simulationbecause the SSN's speed advantage was not sufficient to be able to trap the 


Kilo within the area of diminishing squares. !° 


Table 3 shows the average results at each trial point involving the different search 
paths used, including the random search. With the same input parameters and limitations 
as the random active sonar search, the lawnmower path performance was on par with the 
random search path. This suggests that the Kilo's speed is sufficient to reduce the 


performance of the lawnmower search to be no better than a random search. The spiral- 


10 Washburn, Search and Detection, 1-11. 
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out search path, however, performed better than a random search at every data point. 


Table 4 shows an overall reduction in detection time of 20%. 


Avg Time To 
Avg Time to Init. Detect Avg Time to 
Ping Init. Detect Using Init. Detect 
Det Interval Using Random Lawnmower Using Spiral 
Prob (hours) Search Search Out 


141.19825 
64.4391 

117.37425 
3.32905 
48.61055 
63.48285 


72.67396985 51.6569 49.2787 
87.25525 98.26395 


2.68425 
47.8445 
29.1944 
417.1178 
57.1936 
75,79385 
50.1851 
150.4501 


Overall 
110.74396 122.3584588 87.57010294 


Table 3. | Comparison of Search Paths 


95.82585 90.6327 72.2219 
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Avg Time to Init. Detect | Avg Time to Init. Detect 
sing Random Search Using Spiral Out Percent Improvement 
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35.69445 33.326 663% 
67.2388 27.70% 
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95.82585 ; jo 463% 
28.6538 | 
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87.25525 aa aay 
50.9184 J] 604% 
36.0078 J i892% 
H16.83% 

63.06215 ir a 
75.19385 BP 3212% 
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ae ee B3135% 
yo204] 19.79% | 


Table 4. | Comparison of Random and Spiral-Out Search Paths 


501.518 


Using 





It was previously noted that the Kilo's speed is the factor that causes an 
exhaustive search to decline into a random search's performance. This is due to the Kilo's 
ability to enter previously searched areas while the SSN continues conducting its 
assigned search pattern. In a lawnmower search path, the perimeter that the Kilo can 
cross to enter the previously searched area is the length of the operational area (125nm). 
In the spiral-out search path, this perimeter begins at zero and then increases linearly as 
the search is conducted. As a result, it is difficult in the first part of the simulation for the 
Kilo to enter the previously searched area. This yields a more efficient search as well as 


a reduction of initial detection time. 


In order to compare an active sonar search using a spiral-out technique to the 
passive sonar search, simulation runs were performed to chart the time to initial detection 


using active sonar as a percentage of time to a passive sonar search. The simulations 
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used the same submarine speeds and a ratio of active detection range to passive detection 
range. Figure 8 shows the result of the simulation runs, with the following input 


parameters: 
SSN Speed=15 knots 
Kilo Speed=3 knots 
Passive Detection Range=2nm 
Active detection probability=0.75 
Ping Interval=0.25 hours 


Active Range is varied from 2*(passive det. Range) to 


7*(passive det. range) 
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Figure 8. Active vs. Passive Detection Times 


Initially, when active sonar range is twice that of passive sonar range, the passive 
search performs better. This is because active sonar pings occasionally miss the Kilo 
unless the SSN is constantly pinging or active sonar has a detection probability of 1. As 
active sonar range increases, there is a significant reduction of search times that follow an 
exponential trend. When active sonar range is five times that of passive sonar range, the 
active search detects the Kilo in 25% of the time that a passive search would take. 


Though this is not predictive of the real world situation, the results imply that an active 
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search can detect a submarine in a fraction of the time than a passive search takes. Real- 
world exercises and data could determine the breakpoint at which the active sonar search 


will outperform the passive sonar search. 
D PING INTERVAL ANALYSIS 


It has been shown that initial detection time has an exponential relationship with 
ping interval. This means that as the time between pings increases, there is very little 
initial loss in detection time. This effect builds, however, such that eventually initial 
detection time increases substantially. This is shown in Figure 9. SSN speed was set to 
10 knots,to the Kilo’s 3 knots, active detection range is 10nm, and the probability that a 
ping inside the detection range leads to detection is 0.75. 200 simulations were then 


performed on a series of ping intervals that is increased by 0.25 hours. 


Detection Time vs. Ping Interval 
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Figure 9. Detection Time vs. Ping Interval 


Since initial detection time is not only a function of ping interval, the other factors 
may be adjusted to mitigate the exponential increase in detection time versus ping 
interval. This will deteriorate the best case scenario of initial detection time, but due to 
uncertainties in a real life situation it may be appropriate in order to gain the reduction of 


risk associated with less frequent pinging. 
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1, Ping Interval vs. SSN Speed 


The first, and easiest, parameter for the SSN to change during search is its speed. 
Additional simulations to those performed in section D were conducted--once with SSN 
speed being reduced to 5 knots, and another with search speed increasing to 10 knots. All 


other parameters are held constant. Figure 10 shows the results of the simulations. 
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Figure 10. Ping Interval and SSN Speed vs. Detection Time 


With a very short ping interval, a speed of 15 knots provides the best initial 
detection time at 29.46 hours. Over a ping interval range of 1 to 1.75, the initial detection 
time increases from approximately 71 hours to over 520 hours. Alternately, at 10 knots 
the initial detection time does not increase above 200 hours until a 2 hour ping interval is 
reached, and at 5 knots it is not reached until the ping interval is 2.75 hours. If a 
Commanding Officer feels that a threat exists and is willing to sacrifice the best initial 


detection time for a reduction in ping interval, he should reduce SSN search speed. 
2. Ping Interval vs. Detection Probability 


One factor for a Commanding Officer to consider in his decision to use active 


sonar is the probability of detection if the Kilo is within active sonar range. For 
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example, if the SSN operates in a high clutter environment and the detection probability 
is lowered, it should be understood how much detection time will change if ping interval 


is reduced. 


Figure 11 shows the results of simulation runs similar to those performed in 
Section 1. However, instead of varying SSN speed, detection probability is varied from 
0.4 to 0.9. Varying the detection probability does not have the same impact on detection 
time as SSN speed had. As long as sufficient looks are performed on the Kilo as it is 
within detection range, there will be a moderately high probability of detecting the Kilo. 
For example, assume the SSN is able to ping the Kilo twice as it is in detection range. 


The probability the SSN can detect the Kilo in these two looks is: 


Prob(Detected in two looks)= Prob(Detected in the First Look)+Prob(Missed in First 
Look)*Prob(Detected in Second Look) 


As a result, at a detection probability of 0.4, the overall chance of detecting the 
Kilo in two looks is 64%. At a detection probability of 0.7 and 0.9, this chance increases 
to 91% and 99% respectively. The ping interval in these simulations are generally short 
enough such that multiple looks at the Kilo in detection range is obtained, and detection 


probability does not have as large an affect on detection time as ping interval is varied. 
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Figure 11. Ping Interval and Detect Prob. Vs. Detection Time 
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3. Ping Interval vs. Active Range 


The last factor known to the Commanding Officer is the estimated active 
detection range. Comparably to detection probability, if the active range is considered 
low and the Commanding Officer feels it tactically necessary to reduce ping interval, 
there needs to be an understanding of how this will affect active detection time. Figure 


12 shows the effect of changing ping interval as active detection range changes from 5nm 
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Figure 12. Ping Interval and Active Range vs. Detection Time 


Active range is extremely significant in exacerbating the effects of the exponential 
distribution on detection time. If active range is significant, it takes a drastic delay in ping 
interval to see the exponential effect. Conversely, a relatively short detection range will 
result in a very rapid increase in detection time if ping interval is delayed. Therefore, active 
range is the primary factor the Commanding Officer should consider if he decides to employ 


active sonar. 


Statistically, the model derived in Sections A and B are not a good fit, as they only 
explain approximately 20% of the variability in the model. Further research, including more 
scenarios and possibly additional factors would be needed to attempt to explain more of 
variability in the scenario. The graphs do strongly suggest and exponential relationship of 
the factors in regards to initial detection time, and this relationship still should be included in 


any decision to use active sonar when detecting diesel submarines. 


32 


Vv. CONCLUSIONS 


The following synopsis provides the outcome of the simulations performed in this 


thesis: 


e The performance graphs of the simulations do strong suggest an 
exponential pattern, implying that the time to initially detect a target is an 
exponential function based on detection range, speed of the participants, 
and in the case of active sonar ping, interval and detection probability. As 
a result, as the passive detection range of diesel submarines continues to 
decrease, the average detection time will increase exponentially. There is 
a point where a passive prosecution is no longer feasible, and either other 
assets need to be employed, or an active search could be performed to 


maintain a reasonable detection time. 


e The speed of a submerged diesel submarine is limited; however, it is 
sufficient to make many search patterns no more effective than a random 
search if the diesel can enter previously searched areas. A "spiral-out" 
path can provide an additional decrease in detection time, as it provides an 
initially small perimeter that the diesel submarine can cross in order to 
enter this previously searched area. The simulations suggest that an active 
search can yield detection times at a fraction of passive searches, but real 
world exercises are needed to determine the breakpoint at which this 


Occurs. 


e If the Commanding Officer deems it necessary to reduce ping interval to 
minimize the threat of being counter-detected, the following factors (listed 
in order of priority) will mitigate the effect of the exponential increase in 


detection time: 


1. A long active detection range 
2. Traveling at slower search speeds 
2: A moderate to high detection probability 
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To advance the study of active sonar employment, further research is needed to 
identify additional tactics to optimize its performance. Such studies include identifying a 
set of possible actions once detection occurs. For example, if engagement with the diesel 
submarine is desired, there must be a determination whether to rapidly engage the diesel 
or to continue tracking and obtain better tactical positioning. A proper combination of 


active and passive sonar tracking can also be explored. 


The regression model performed in this thesis only explains approximately 20% 
of the variability that appears in the simulations. As a result, further research is needed to 
attempt to better explain this variability. These may include additional scenarios and 
additional factors. Future work in this field should include an attempt to build a better 
statistical model that better explains the variability shown. Additional scenarios could 
include factors that were not considered in this thesis, and may also contain situations 
that would not occur in a real world encounter. An incremental increase in the number of 
replications could identify whether the randomness of the model or the situation being 


modeled overcomes the factors due to its stochastic nature. 


The modern diesel submarine provides significant challenges to the U.S. naval 
fleet. Not only does it pose a significant threat to high value units such as aircraft 
carriers, but the relatively cheap cost of each diesel submarine is such that they may be 
deployed in significant numbers to overwhelm small numbers of high quality units such 
as nuclear submarines. The ability to detect diesel submarines early provides the best 


ability to successfully defend against any threat that may be posed. 
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APPENDIX 


public class ActiveSearchRandom { 


as 
* The class random number generator will produce a different random number 

* at the start of each run. The parameters are the number of replications, 

* SSN speed, Diesel Speed, Active detection range, probability of 

* detection, and ping interval which can be overridden using the arguments. 

* 

* @param args The values of the number of replications run, SSN speed, Diesel Speed, Active 
*Detection Range, Detection Probability, Ping Interval, 

yf 


public static void main(String args[]) { 


double simTime = 0.0; 
double totalSimTime = 0.0; 
int numberReplications = 200; 
double ssnSpeed = 0.0; 
double dieselSpeed = 0.0; 
double activeRange = 0.0; 
double probDetect = .0; 
double pingInterval = .0; 
if (args.length == 6) { 
numberReplications = Integer.parse/nt(args[0]); 
ssnSpeed = Double.parseDouble(args[1]); 
dieselSpeed = Double.parseDouble(args[2]); 
activeRange = Double.parseDouble(args[3]); 
probDetect = Double.parseDouble(args[4]); 
pingInterval = Double.parseDouble(args[5]); 
} 
SensorTargetMediatorFactory.addMediator(ActiveSonarSensor2.class, 
UniformLinearMover.class, ActiveSonarMediator2.class); 


RandomVariate randPosit = RandomVariateFactory.getinstance( 
"Uniform", new Object[] { new Double(0.0), 
new Double(125.0) }); 
Point2D.Double ssnStart = new Point2D.Double(randPosit 
-generate(), randPosit.generate()); 
Point2D.Double dieselStart = new Point2D.Double(randPosit 
-generate(), randPosit.generate()); 


// Create the diesel and the SSN as uniform linear movers. 

Mover badGuy = new UniformLinearMover("Diesel", dieselStart, 
dieselSpeed); 

Mover goodGuy = new UniformLinearMover("SSN", ssnStart, ssnSpeed); 


// Attach the passive sonar to the SSN 
Sensor ssnActive = new ActiveSonarSensor2(goodGuy, activeRange, 


25 


pingInterval, probDetect); 


// Register the movers to the Referee 

SensorTargetReferee activeRef = new SensorTargetReferee(); 
activeRef.register(badGuy); 

activeRef.register(ssnActive); 


goodGuy.moveTo(new Point2D.Double(randPosit.generate(), randPosit 
.generate())); 

badGuy.moveTo(new Point2D.Double(randPosit.generate(), randPosit 
.generate())); 

// Have the Diesel and SSN to move randomly 


RandomVariate[] rv = new RandomVariate[] { 
RandomVariateFactory.getInstance(""Uniform", new Object[] 


0.0, 
125.0 }), 
RandomVariateFactory.getInstance(""Uniform", new Object[] 


0.0, 
125.0 })}; 


RandomLocationMoverManager goodGuyMM = new 
RandomLocationMoverManager( 
goodGuy, rv); 
RandomLocationMoverManager badGuyMM = new 
RandomLocationMoverManager( 
badGuy, rv); 


goodGuyMM.setStartOnReset(true); 

badGuyMM.setStartOnReset(true); 

SandboxFrame frame = new SandboxFrame("Passive Sensor Platform"); 

frame.setSize(1000, 1000); 

((PingThread) frame.getControlPanel().getController()) 
.setMillisPerSimtime(50); 


frame.getSandbox().setOrigin(new Point2D.Double(50, 150)); 
frame.addMover(badGuy, Color.red); 
frame.addMover(goodGuy, Color.blue); 


frame.addSensor(ssnActive, Color.orange); 
PropertyChangeFrame pcf = new PropertyChangeFrame(); 


ssnActive.addPropertyChangeListener("detection", pcf); 
ssnActive.addPropertyChangeListener("undetection", pcf); 


frame.setLocation(1, 1); 

frame.setVisible(true); 

pef.setLocation(frame.getLocationOnScreen().x + frame.getWidthQ), 
frame.getLocationOnScreen().y); 

pcf.set Visible(true); 
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for (int i=0;i<numberReplications;i++) { 
Schedule.resert(); 

Schedule.setVerbose(false); 
Schedule.stopOnEvent(1, "Detection", Contact.class); 


Schedule.startSimulation(); 
simTime = Schedule.getSimTime(); 


DecimalFormat form = new DecimalFormat("##.##"); 
System.out.println(form.format(simTime)); 


} 
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public class PassiveSearchRandom { 


[re 


* The class random number generator will produce a different random number 
* at the start of each run. The parameters are the number of replications, 

* SSN speed, Diesel Speed, Passive detection range, and probability of 

* detection, and can be overriden using the arguments. 


Ba 


* @param args The values of the number of replications, SSN 
*Speed, Diesel Speed, and Passive sonar detection range. 


yh 


public static void main(String args[]) { 


double simTime = 0.0; 

int numberReplications = 200; 

double ssnSpeed =15.0; 

double dieselSpeed = 3.0; 

double passiveRange = 2.; 

if (args.length == 4) { 
numberReplications = Integer.parseInt(args[0]); 
ssnSpeed = Double.parseDouble(args[1]); 
dieselSpeed = Double.parseDouble(args[2]); 
passiveRange = Double.parseDouble(args[3]); 


SensorTargetMediatorFactory.addMediator(CookieCutterSensor.class, 
UniformLinearMover.class, CookieCutterMediator.class); 
// Establishing a RandomVariate position for the Diesel and SSN to 
// be generated. 
RandomVariate randPosit = RandomVariateFactory.getinstance( 
"Uniform", new Object[] { new Double(0.0), 
new Double(125.0) }); 


Point2D.Double ssnStart = new Point2D.Double(0.0, 0.0); 
Point2D.Double dieselStart = new Point2D.Double(randPosit 
-generate(), randPosit.generate()); 
// Create the diesel and the SSN as uniform linear movers. 
Mover badGuy = new UniformLinearMover("Diesel", dieselStart, 
dieselSpeed); 
Mover goodGuy = new UniformLinearMover("SSN", ssnStart, ssnSpeed); 


// Attach the passive sonar to the SSN 
Sensor ssnPassive = new CookieCutterSensor(passiveRange, goodGuy); 


// Register the movers to the Referee 

SensorTargetReferee passiveRef = new SensorTargetReferee(); 
passiveRef.register(badGuy); 

passiveRef.register(ssnPassive); 


goodGuy.moveTo(new Point2D.Double(0.0, 0.0)); 
badGuy.moveTo(new Point2D.Double(randPosit.generate(), randPosit 
.generate())); 


// Have the Diesel and SSN to move randomly 
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RandomVariate[] rv = new RandomVariate[] { 
RandomVariateFactory.getInstance(""Uniform", new Object[] 


0.0, 
125.0 }), 
RandomVariateFactory.getInstance(""Uniform", new Object[] 


0.0, 
125.0 })}; 


// RandomLocationMoverManager goodGuyMM=new 
// RandomLocationMoverManager(goodGuy, rv); 
PatrolMoverManager goodGuyMM = new PatrolMoverManager(goodGuy, 
pathlawn); 
RandomLocationMoverManager badGuyMM = new 
RandomLocationMoverManager( 
badGuy, rv); 


goodGuyMM.setStartOnReset(true); 
badGuyMM.setStartOnReset(true); 


SandboxFrame frame = new SandboxFrame("Passive Sensor Platform"); 

frame.setSize(150, 150); 

((PingThread) frame.getControlPanel().getController()) 
.setMillisPerSimtime(50); 


frame.getSandbox().setOrigin(new Point2D.Double(0, 160)); 
frame.addMover(badGuy, Color.red); 
frame.addMover(goodGuy, Color.blue); 


frame.addSensor(ssnPassive, Color.orange); 
PropertyChangeFrame pcf = new PropertyChangeFrame(); 


ssnPassive.addPropertyChangeListener("detection", pcf); 
ssnPassive.addPropertyChangeListener("undetection", pcf); 


frame.setLocation(1, 1); 

frame.setVisible(true); 

pcf.setLocation(frame.getLocationOnScreen().x + frame.getWidth(), 
frame.getLocationOnScreen().y); 

pcf.set Visible(true); 


for (int i=0;i<numberReplications; i++) { 
Schedule.resetQ); 

Schedule.setVerbose(false); 
Schedule.stopOnEvent(1, "Detection", Contact.class); 


Schedule.startSimulation(); 

simTime = Schedule.getSimTime(); 

DecimalFormat form = new DecimalFormat("##.##"); 
System.out.printIn(form.format(simTime)); 


} 
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RSquare 0.205821 
RSquare Adj 0.205119 
Root Mean Square Error 752.0049 
Mean of Response 643.015 
Observations (or Sum Vgts) 3400 
| Analysis of Variance 
Sum of 

Source DF Squares Mean Square F Ratio 
Model 3 497713105 165904368 293.3705 
Error 3396 1920476672 $65511.39 Prob>F 
Cc. Total 3399 2418189777 <.0001* 
| Lack Of Fit 

Sum of F Ratio 
Source DF Squares Mean Square 18.8109 
Lack Of Fit 13 129464409 9958801 Prob>F 
Pure Error 3383 1791012262 $29415 <.0001* 
Total Error 3396 1920476672 Max RSq 

0.2594 

| Parameter Estimates 
Term Estimate Std Error tRatio Prob>jt| 
intercept 2018.1492 66.0328 30.56 <.0001* 
SSN Spd (hours) 64.0973 4.255024 -15.06 <.0001* 
Kilo Spd (hours) 9.7152394 13.5023 0.72 0.4719 
Pass Rng (nm) -436.2267 17.33581 -25.16 <.0001* 
| Sorted Parameter Estimates 
Term Estimate StdeError t Ratio Prob>jt| 
Pass Rng (nm) 436.2267 17.33581 -25.16 <.0001* 
SSN Spd (hours) 64.0973 4.255024 -15.06 ES <.0001* 
Kilo Spd (hours) 9.7152394 13.5023 0.72 0.4719 
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Figure 13. Linear Passive Sonar Regression 


40 





| Response LN (Avg Time Init Det) 
| Actual by Predicted Plot 











Init Det) Actual 


2 
£ 
— 

> 
= 
=z 
=< 


3.0 40 50 60 70 80 9.0 
LN (Avg Time Init Det) Predicted 
P<.0001 RSq=0.23 RMSE=0.9718 








| Summary of Fit 





























RSquare 0.232445 
RSquare Adj 0.231767 
Root Mean Square Error 0.971785 
Mean of Response 5.87993 
Observations (or Sum Vvgts)} 3400 
| Analysis of Variance 
Sum of 

Source DF Squares Mean Square F Ratio 
Model = | 971.2218 323.741 342.8128 
Error 3396 3207.0658 0.944 Prob>F 
Cc. Total 3399 4178.2876 <.0001* 
| Lack Of Fit 

Sum of F Ratio 
Source DF Squares Mean Square 4.17859 
Lack Of Fit 13 50.6505 3.89612 Prob>F 
Pure Error 3383 3156.4153 0.93302 <.0001* 
Total Error 3396 3207.0658 Max RSq 

0.2446 

| Parameter Estimates 
Term Estimate StdError tRatio Prob>/jt|/ 
intercept 7.9734544 0.085331 93.50 0.0000* 
SSN Spd (hours) -0.100753 0.005499 -18.32 <.0001* 
Kilo Spd (hours) -0.025637 0.017448 -1.47 0.1418 
Pas= Rng (nm) -0.567381 0.022402 -25.33 <.0001* 
| Sorted Parameter Estimates 
Term Estimate StdError t Ratio Prob>|t| 
Pass Rng (nm) -0.567381 0.022402 -25.33 <.0001* 
SSN Spd (hours) -0.100753 0.005499 -13.32 <.0001* 
Kilo Spd (hours) -0.025637 0.017448 -1.47 Eee 0.1418 
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Correlations 

LN (Avg Time Init Det) SSN Spd (hours) Kilo Spd (hours) Pass Rng (nm) 
LN (Avg Time Init Det) 1.0000 -0.2782 -0.1193 -0.3935 
SSN Spd (hours) -0.2782 1.0000 0.0688 0.0012 
Kilo Spd (hours) -0.1193 0.0688 1.0000 0.2000 
Pass Rng (nm) -0.3935 0.0012 0.2000 1.0000 
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Figure 15. Multivariate Analysis in Log-Linear Passive Model 
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RSquare 0.168596 
RSquare Adj 0.167371 
Root Mean Square Error 194.8522 
Mean of Response 124.9744 
Observations (or Sum Wots) 3400 
| Analysis of Variance 
Sum of 

Source DF Squares Mean Square F Ratio 
Model 5 26131046 5226209 137.6500 
Error 3394 128861261 37967 Prob>F 
C. Total 3399 154992308 <.0001* 
| Lack Of Fit 

Sum of F Ratio 
Source DF Squares Mean Square 65.3386 
Lack Of Fit 11 22579718 2052702 Prob>F 
Pure Error 3383 106281544 31416 <.0001* 
Total Error 3394 128861261 Max RSq 

0.3143 

| Parameter Estimates } 
Term Estimate StdError tRatio Prob>|t| 
Intercept 319.10836 23.92878 13.34 <.0001* 
SSN Speed (knots) -4.827698 1.113139 -4.34 <.0001* 
Kilo Speed (knots) 18.306534 3.4785 §.26 <.0001* 
Active Range (nm) -22.78355 1.092225 -20.86 <.0001* 
Det Prob -83.71544 15.57372 -§.70 <.0001* 
Ping Intw 149.688 11.06756 13.52 <.0001* 
| Sorted Parameter Estimates 
Term Estimate StdError t Ratio 
Active Range (nm) -22.78355 1.092225 -20.86 
Ping Intw 149.688 11.06756 13.52 
Det Prob -88.71544 15.57372 -5.70 
Kilo Speed (knots) 18.306534 3.4785 5.26 
SSN Speed (knots) -4.827698 1.113139 -4.34 
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RSquare 0.2005 
RSquare Adj 0.199322 
Root Mean Square Error 1.328495 
Mean of Response 3.91214 
Observations (or Sum Wgts) 3400 
| Analysis of Variance 
Sum of 

Source DF Squares Mean Square F Ratio 
Model 5 1502.198S 300.440 170.2305 
Error 3394 $990.0675 1.765 Prob>F 
C. Total 3399 7492.2660 <.0001* 
| Lack Of Fit 

Sum of F Ratio 
Source DF Squares Mean Square 23.7739 
Lack Of Fit 11 429.9201 39.0836 Prob>F 
Pure Error 3383 $560.1473 1.6436 <.0001* 
Total Error 3394 5990.0675 Max RSq 

0.2579 

| Parameter Estimates 
Term Estimate StdError tRatio Prob>jt] 
intercept 6.078588 0.163146 37.26 <.0001* 
SSN Speed (knots) 0.05499 0.007589 7.25 <.0001* 
Kilo Speed (knots) 0.0474707 0.023716 2.00 0.0454" 
Active Range (nm) -0.173232 0.007447 -23.26 <.0001* 
Det Prob -0.86098S 0.106181 3.11 <.0001* 
Ping Intw 1.0165373 0.075458 13.47 <.0001* 
| Sorted Parameter Estimates } 
Term Estimate StdError t Ratio Prob>|t| 
Active Range (nm) -0.173232 0.007447 -23.26 <.0001* 
Ping Intv 1.0165373 0.075458 13.47 <.0001* 
Det Prob -0.360985 0.106181 8.11 <.0001* 
SSN Speed (knots) -0.05499 0.007589 7.25 <.0001* 
Kilo Speed (knots) 0.0474707 0.023716 2.00 0.0454" 
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Figure 17. Log-Linear Active Sonar Regression 
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_ Multivariate 








Correlations 
SSN Speed (knots) Kilo Speed (knots) Active Range (nm) Det Prob 

SSN Speed (knots) 1.0000 -0.1492 0.0233 0.0270 0.0233 
Kilo Speed (knots) -0.1492 1.0000 0.0470 -0.0267 -0.0948 
Active Range (nm) 0.0233 0.0470 1.0000 -0.0028  -0.0004 
Det Prob 0.0270 -0.0267 -0.0028 1.0000 0.0001 
Ping Intv 0.0233 -0.0948 -0.0004 0.0001 1.0000 
LN (Initial Detection Time) -0.1241 0.0149 0.3585 -0.1274 0.2022 
| Scatterplot Matrix 

Sineesr* 

heteertan ine 
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Figure 18. Multivariate Analysis in Log-Linear Active Model 
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Table 5. Ping Interval and Speed vs. Detection Time 
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89.9396 
136.03525 
222.1954 
Poa 562.2488 
a ee . 

1 : 


ae 95.6014 50.8127 
157.3704 54.1943 


——) a 
Po 108.4174 
a: a p34. 7461 


Table 6. _ Ping Interval and Active Detection Range vs. Detection Time 
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1.25 73.5598 42.4747 


5 

1 . 

5 

2 
5 : 

3 





BIBLIOGRAPHY 


"Briefing on the Virginia Class Submarine Contract." 
http://www.defenselink.mil/transcripts/transcript.aspx ?transcriptid=3224 
(accessed Feb 2, 2008). 


Cioppa, Thomas M. "Efficient Nearly Orthogonal and Space-Filling Experimental 
Designs for High-Dimensional Complex Models." Master's Thesis, Naval 
Postgraduate School, 2002. 


Devore, Jay L. "Regression with Transformed Variables." Probability and Statistics for 
Engineering and the Sciences. 6th ed., 565-568. Belmont, Ca: Brooks/Cole- 
Thomson Learning, 2004. 


"Kilo Class Submarine - People's Liberation Army Navy." 
http://www.globalsecurity.org/military/world/china/kilo.htm (accessed Feb 2, 
2008). 


"Song Class - People's Liberation Army Navy." 
http://www. globalsecurity.org/military/world/china/song.htm (accessed Feb 2, 
2008). 


Saunders, Stephen. "U.S. Submarines." Jane's Fighting Ships 2007-2008. Surrey, U.K.: 
Jane's Information Group Limited, 2007. 


"United States Navy Fact File: Attack Submarines-SSN." 
http://www.navy.mil/navydata/fact_display.asp?cid=4100&tid=100&ct=4 
(accessed February 2, 2008). 


Washburn, Alan R. "Models Based on Detection Rate" Search and Detection. 4th ed., 


1.11-2.6. Linthicum, Md: Institute for Operations Research and the Management 
Sciences, 2002. 


47 


THIS PAGE INTENTIONALLY LEFT BLANK 


48 


10. 


11. 


INITIAL DISTRIBUTION LIST 


Defense Technical Information Center 
Ft. Belvoir, Virginia 


Dudley Knox Library 
Naval Postgraduate School 
Monterey, California 


Arnold Buss 
Naval Postgraduate School 
Monterey, California 


Thomas Lucas 
Naval Postgraduate School 
Monterey, California 


James Eagle 
Naval Postgraduate School 
Monterey, California 


RADM Raymond Jones (ret) 
Naval Postgraduate School 
Monterey, California 


CAPT Doug Otte 
Naval Postgraduate School 
Monterey, California 


CAPT Wayne Hughes 
Naval Postgraduate School 
Monterey, California 


Jeff Kline 
Naval Postgraduate School 
Monterey, California 


Pierre Corriveau 
Naval Postgraduate School 
Monterey, California 


Peter Lorenz 
Naval Warfare Development Center 
Newport, Rhode Island 


49 


ED: 


13. 


14. 


CAPT Lee Rosetti 
Naval Warfare Development Center 
Newport, Rhode Island 


Jim Miller 
Submarine Development Squadron 12 
New London, Connecticut 


CAPT Mike McLaughlin 


COMSUBFOR 
Norfolk, Virginia 


50 


